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INVISCID ANALYSIS OF JET INJECTION BETWEEN TWO MOVING STREAMS 
by Marvin E. Goldstein and Robert Siegel 
Lewis Research Center 

SUMMARY 

An analytical method was developed for determining the flow characteristics of a 
two-dimensional jet injected between two moving streams. The jet is flowing out of a 
channel and is turned as it enters into the main stream flow. As a result, the flow field 
consists of three regions: the jet, and the main stream on either side of it. In the solu- 
tion, these regions must be coupled along the two dividing streamline boundaries. The 
velocity variations along these streamlines produce static pressure variations along 
them, which introduces a complexity not considered in the usual free streamline analysis 
of jets. Both external streams have the same total pressure, which is different from 
that in the jet. The solution is obtained for conditions when the total pressure in the jet 
does not differ from the free -stream value by a large amount compared with the free- 
stream dynamic head. Results are obtained for the shape of the jet boimding streamlines 
as a function of the injection geometry and as a function of the total pressure in the jet 
relative to that in the free stream. 


INTRODUCTION 

There are various situations where a stream of fluid is injected into a moving stream 
of a similar fluid (fig. 1). Examples are the discharge of water into a river, gas blown 
from an exhaust into the wind, and the injection of gaseous fuel into a stream of oxidant. 
If, for purposes of analysis, the configuration is approximated as two dimensional, the 
jet region is bounded on both sides by the outer stream flow. Any turning in the flow pro- 
duces a variable pressure along both jet botmdaries, a feature not taken into accoimt in 
classical free-jet analyses, where the external flow is usually at rest and hence at uni- 
form static pressure. In reference 1 an inviscid two-dimensional jet was considered 
exiting from a slot and interacting with a crossflow aloi^ one of the jet boundaries. The 
present report is a further development of the anal 3 rtical method to include a matching of 
a jet with an external flow aloi^ both jet boundaries. As a result of the inviscid 



Outer 

stream 



(a) Flow turned and contracted. 



^^^^Outer 

(b) Flow turned and expanded. 

Figure 1. - Flow configuration. 

assumption the flow patterns are expected to be valid near the jet outlet, where mixing 
effects are small. 

The analysis also accounts for the total pressure in the jet being different from the 
total pressure in the outer stream. The difference between these total pressures must, 
however, not be too large compared with the dynamic pressure head of the undisturbed 
external stream. The problem is solved by expanding the solutions asymptotically in this 
normalized difference in total pressures. 

The zeroth-order solution in the asymptotic expansion corresponds to equal total 
pressures in jet and stream, and the corresponding solution is obtained by a conventional 
conformal mapping method. When the total pressures are different, the jet boundaries 
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will shift away from the zeroth-order streamlines. In the solution method the conditions 
alor^ these boimdaries are transformed to the known zeroth-order boundaries. 

The flow configurations studied allow for tvirning and either contractile or expanding 
the injected stream as it enters the main stream flow, as shown in figures 1(a) and (b). 
Results are given for the shape of the jet free streamlines for these configurations for 
various values of such geometric parameters as the amount of turning. 


ANALYSIS 

Formulation and Boundary Conditions 

It will be assumed that the flow is inviscid, incompressible, and irrotational. The 
notation describing the jet configxiration is illustrated in figure 2(a), which corresponds 




(b) Dimensionless physical plane (z-plane). 
Figure 2. - Flow geometry and nomenclature. 
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to figure 1(a). For the geometry in figure 1(b) the plate angle a would be negative. 

The analysis is limited to the condition where the difference between the total pressure in 
the jet Pj and the total pressure in the main stream is not too large, so that 
[e [ « 1 where 


P. - P 
= _J 

1 v2 


(1) 


and p is the density of the fluid, and is the velocity of the main stream at infinity. 
(All symbols are defined in the appendix. ) 

Let b be a convenient reference length which will be specified in the course of the 
analysis. The X- and Y -components of the velocity, U and V, respectively, are made 
dimensionless by V^; and the stream function ^ and the velocity potential €> are made 
dimensionless by V^b. Thus, the dimensionless quantities u, v, i//, and <p are defined 
by 


u 



V 





V^b 


cp = 


4> 


OO 


The dimensionless complex conjugate velocity ^ and the dimensionless complex poten- 
tial w are defined, as usual, by 


? = u - iv 


and 


w = + ii// 

With all ler^hs made dimensionless by b (i. e. , x = X/b, y = Y/b, I = L/b, 6 = A/b, 
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h = H/b, and r = R/b), the flow configuration is shown in the dimensionless physical 
plane (with the complex variable z defined by z = x + iy) in figure 2(b). 

The stream of fluid issuing from the orifice formed by the walls 345 and 56 
(fig. 2(b)) meets the main stream at the points 3 and 6 and forms the common stream- 
lines which are denoted by and in figure 2. Points on the streamline are 
denoted by z^^^ = x^^^ + iy^^^, and points on the streamline are denoted by 

z(2) = ^(2) ^ iy(2) 

In order to satisfy the requirement that there be no discontinuities in static pressure 
an 5 Twhere within the flow field, it is necessary (as shown in ref. 1) to allow the velocity 
to be discontinuous across and For this reason the streamlines and 

will be called slip lines. The region within the jet and the orifice is denoted by D^, and 
the remaining region of the flow (i. e. , the main stream) by D^. Since the velocity (and, 
as a consequence, the velocity potential) is discontinuous across and it is 

convenient to use a superscript I to denote the flow quantities inside the jet (i. e. , in D^) 
and a superscript II to denote those in the main stream. Thus, 


and 



r^^(z) 

for 

z C 

C(z) 

1 

,?°(z) 

for 

zCD^I 

i 

r T 
w^(z) 

for 

zC 

w(z) 

w°(z) 

for 

zC 


Then and w^ are hoi omorphic in the interior of and and w^ are holomor- 
phic in the interior of D^. 

A repetition of the argument given in reference 1 shows that Bernoulli’s equation 
implies 




2 _ Pj - P, 


1 ,j2 


= € 


for every point z^^^ of and 


. (5n(,(2))p , 


( 2 ) 


(3) 
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for every point of Since and S^^^ are both common streamlines, it is 

clear that 

9m w'(J«) = 9m. w“(z»>) = constant 
9m w'(z<2>) = 9m w“(z®) = constant 

Moreover, the arbitrariness in the level of w can be partially removed by choosing 
one of these constants to be zero (ref. 2). Hence, we can put 

9m = 9m w^(z^^^) = 0 (4) 

(5) 

12) 

where p' is a real constant. 

The remaining arbitrariness in w can be removed by choosing (at the end point of 
the flap in fig. 2(b)) 


w^(Ze"^’^^) = w°(Ze“^’^^) = 0 (6) 

The condition imposed on the velocity at infinity is (in view of the manner of nondimen - 
sionalization) 


^^(z) — 1 for z ^ point 1 and z — point 7 (7) 

The remaining boundary conditions are that the normal component of the velocity 
vanish on the solid boundaries. These conditions are sufficient to completely determine 
the solution. They are summarized here for convenient reference: 
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|?\z)| -|?«(z)| =e 

9m w^(z) = 9m w^(z) = 0^ 


for z C S 


(1) 


?^Z) 


-k"(z)l =e 


for z C S 


( 2 ) 


9m w^(z) = 9m w^(z) = >- 

9m ?^(z) =0 for z C 4? and z C 56 

9m ?^(z) = 0 for z C 12 and z C 67 

arg ?^(z) = ttX. for z C 34 

arg ?^^(z) = 7T\ for z C 23 


(8) 


Asymptotic Expansions 

For small values of e the functions w\ and can be expanded in an 

asymptotic power series in e. In view of the fact that the shape of each slip lines de- 
pends on e, these expansions imply that the coordinates of as well 

as the asymptotic jet width 6, must also be expanded in powers of e. Hence, 


+ • • • 

+ . . . 

= Wq + + . . . 

= Wq + + . . . 

= z^^^ + + . . . 

6 = + . . . 

^( 2 ) = ^ 2 ) ^ ,^( 2 ) ^ _ 


(9) 
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As pointed out in reference 1 the expansion of and does not imply that the 
complex variable z itself is being expanded. It is shown in reference 1 that the coeffi- 
cients in the expansions of and are related to those in the expansions of and 
w^, respectively, by 


( 10 ) 


The reason for omitting the superscript I or n in the zeroth -order terms of the 
first four expansions is that (as will be shown subsequently) the zeroth-order solutions 
are not discontinuous across the curves and Hence, there is a single function 

which is holomorphic in the entire flow field (of course, the same is true of Wq). 

The reference length b will now be chosen in such a way that 6 q = 1. Thus, b is 
the zeroth-order asymptotic thickness (at «>) of the jet. This is denoted symbolically by 
putting b = Aq. The second last expansion (9) is then 

6 = 1 + e5j + O(e^) . . . (11) 

Zeroth-Order Solution 

Boundary conditions for zeroth-order solution . - When the expansions (eq. (9)) are 
substituted into the boundary conditions (eqs. (7) and (8)) and only the zeroth-order terms 
are retained, the following bo\andary conditions for the zeroth-order solution are obtained. 
First, the first two boxmdary conditions in equation (8) show, as has already been antic- 
ipated, that the zeroth-order solution must be continuous across the slip lines and, hence, 
is characterized by functions which are holomorphic everywhere within the flow field. 

The remaining conditions show that 
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I 


II 

O 


(12) 

U Wq(zW) = oj 

Wo(z[,2)) = ^(2) 

(13) 

^ €45, zC 56, z C 12, and z C 67 

(14) 

arg ?q(z) = 7tX for z C 23 and z C 34 

(15) 

?q(z) — 1 for z — point 1 and z —point 7 

(16) 

The conditions (12) and (13) merely serve to identify the coordinates 

z^^^ and of 


the zeroth-order streamlines emanating from points 3 and 6 in figure 2(b) and to indicate 
that the streamline emanating from point 3 is to be taken as the zero streamline. 

Now the change in the stream function across the jet must be equal to the volume flow 

rate through the jet. Hence, if denotes the dimensionless zeroth-order volume flow 

^ ( 2 ) 

through the jet, and since equation (13) shows that -\p^ ' is the zeroth-order change in 
the stream function across the jet, it is clear from the definition of the stream function 
that Qq = 

However, the boxmdary condition (16) shows that far downstream in the jet the zeroth- 
order velocity goes to 1, and equation (11) shows that the asymptotic thickness of the jet 
must also be 1. It follows from these remarks that Qq = 1. Hence, the last expan- 
sion (9) becomes 

= -1 + + . . . (17) 

Representation in potential and hodograph planes. - Now the boundary value problem 
posed by the boundary conditions (14) to (16) is a simple free streamline problem which 
can be readily solved by the Helmholtz -Kir chhoff technique. The procedure for obtaining 
the solution is (ref. 2) to draw the region of flow in the hodograph plane and in the com- 
plex potential plane, and then to find the appropriate mapping of these two planes into 
some convenient intermediate plane (say, the T -plane). The shapes of these regions can 
readily be deduced from the boundary conditions (14) to (16), and they are shown in fig- 
ures 3 and 4. (We have put Wq = + ii/^Q in fig, 3, ) The corresponding points in the 

various planes are designated by the same numbers. The zeroth-order "slip lines" are 
shown dashed in these figures since they do not correspond to lines of discontinuity and 
therefore can be ignored as far as obtaining the zeroth-order solution is concerned. 
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Figure 3. - Zeroth-order complex potential plane (WQ-plane). 



Figure 4. - Zeroth-order hodograph (^-plane). 


Intermediate T -plane. - The intermediate T -plane is chosen in such a way that the 
region of flow maps into the upper half plane in the manner indicated in figure 5. We 
shall denote the real and imaginary parts of the variable T by ^ and -q, respectively. 
The region of the T -plane into which the zeroth-order flow field interior to the jet maps 
is denoted by 9^, and the region of the T -plane into which the zeroth-order main stream 
maps is denoted by 9^^. The dividing lines between these two regions corresponding to 
the slip lines and (which for convenience are called the zeroth-order slip lines 
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even though no slip occurs in the zeroth-order solution) are denoted by and 
respectively. The figure is drawn for the case where > y^. The situation where 
^3 ^ ^6 3-lso occur. However, we show only the former case. In the latter case the 
two zeroth-order slip lines will approach each other asymptotically at the left of the imag- 
inary axis. Otherwise the figures will look the same. 

Mappii^ of potential and hodograph planes into intermediate T -plane. - A simple 
application of the Schwartz -Christoffel transformation (ref. 3) shows that the mapping 
which properly transforms the Wq - plane into the upper half T-plane in the manner indi- 
cated in the figures is defined by 


dw. 


dT 


1 <T ^ rgKT - rg) 
iryjTg T 


for ?7 > 0 


or, by performing the indicated integration. 


^0 ■- 

7T 


InJL 

L ^3 



(^6 - ^ 3 ^ 
^^ 3^6 


(T + rg) 


(18) 


(19) 
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Figure 6. - Intermediate plane (Q -plane). 


The mapping of the zeroth-order hodograph plane (fig. 4) into the upper half T-plane 
is best accomplished by introducing the intermediate ^ -plane (with 12 = a + ir) shown in 
figure 6. The mapping that transforms the plane into the upper half 12 -plane in the 
manner indicated in the figxires is 

Cq = for t > 0 (20) 

On the other hand, the linear fractional transform 



transforms the upper half T-plane into the upper half -plane in the manner indicated in 
the figures. The relation between the y’s and /3*s is 
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1 -yg 


^^5 = ^4 > 



^6 + ^4 

re + 1 


( 22 ) 


Equations (20) and (21) can now be combined to show that the mapping of the upper 
half T -plane into the -plane is 


fT + 




74 ' 


vT + 1 


for 77 > 0 


(23) 


Relation of mapped coordinates to physica l dimensions. - It is necessary to relate 
the y's, or equivalently the iS's, to the physical parameters involved in the problem. 
This is most conveniently done by working in the -plane. In order to do this, however, 
it is first necessary to obtain dwg/dfi as a function of To this end notice that equa- 
tion (21) can be solved for T to obtain 


T = 


74-^2 


- 1 


Hence, 


(J2 - 1)^ 

Evidently, 

dJ2 dT dn 


(24) 


(25) 


(26) 
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Substituting in equations (18), (22), (24), and (25) shows that 

dwQ (1 - - y3)(l + y^) (n + ^g)(n - /3g) 

^’^2?'% (fi - - 1)^ 

But equations (22) can be solved for the y's in terms of the /3’s to obtain 


/3o + /3c 


^3 


1 + 




> 


^6-^5 


1 - i^e 


Substituting these into equation (27) gives 

dwQ (/Bg - l)^ (fl + - 0 q ) 

dfi TrOg + /3g)(/3g - /3g) _ ^^3 


But the first equation (10) shows that 


df2 dWg dJ2 dS2 

Hence, it follows from equations (20) and (29) that 



($7 + ^o)(J2 “ /3fi) 

^ ^ do, 

(n - /3g)(fi - 1)^ 


( 27 ) 


(28) 


(29) 


(30) 


The integral can now be used to relate the various lengths in the physical plane 
(fig. 2) to the y’s and/or the /3’s. 

Thus, it can be seen from figure 2(b) that the length I is given by the absolute value 
of the integral carried out between points 3 and 4. Or upon deforming the contour to lie 
along the real axis in the n -plane. 
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I 


Jt (^3 + ^5)^0 - 



. (a + /3^)(a- - /3«) 

-J " 2_dCT 

(<r)^ (cr - ^g)(a - 1)^ 


Upon charging the variable of integration from a to -a this becomes 


I 


(' - Hf 

+ ^5)(/3g - ^5) 



. (/3„ - CT)(a + j3«) 

2 ^ L dcr 

(a + l)^(a + /3g) 


(31) 


On the other hand, in order that the points 2 and 4 correspond to the same point in 
the physical plane, we see from equation (30) that we must have 


0 = 



(ct + i 3 „)(a - / 3 „) 

— ^ da 

(a - /3g)(a - 1)^ 


Expanding the integrand in partial fractions gives 



By changing the variable of integration from a to in the first integral, and from 

a to -a in the remaining integrals, each of these integrals becomes a beta -function. 
We, therefore, obtain 
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Using the properties of the gamma -function to simplify we get, upon dividing through by 
r(X) and r(l - X) and rearranging 

(% - + < 33 ) 1 - -i - X(/3j - 1) = fej - A2(fl3 + Dtfg - 1) - x (^5 - l)3 (32) 

L J ' 

Finally, it can be seen from figure 2(b) that carrying out the integral (30) from 
point 4 to point 6 yields -(r + ih). Hence, upon deforming the contour of integration to 



Figure 7. - Path of integration in Q -plane. 
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the path shown in figure 7 we get in the limit as Cq — 0 




-(r + ih) = 


(' - ^ 5 )= 






^ (/3g + CT ){^Q - cr ) 

CT^(a -/3j(l-a)3 


da 


+ lim 
€0-^ T7(/33 


(n - 


- /Sg) 




- - ir 


where 'f. denotes that the Cauchy principle value is to be taken. Hence, equating 
real and imaginary parts gives 


(' - 




1 (/3q + 0 ^ )(cr - /3g) 

I J__3 ^da 

a^ (a - I3^){1 - (J)^ 


and 


-ih = lim 


(%- 


€q-0 + /3g)(^g - 



+ /3g)(^2 - ( 3 q ) 

- /3g)(S2 - 1)^ 


dJ2 


(33) 


lim 



-L- lim 

7t4 ^0"^ 


In (a - 13^) 




iO 




ITT 


lim [In ^ - In ^ - iTT ] = - -i- 
7T/3J "0“^ 0 I 
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Hence, 


h 



(34) 


Thus, equations (31), (32), and (33) serve to determine /Sg, /3g, and (3^ in terms of the 
lengths I and r; and equation (34) gives the ratio of nozzle width to zeroth-order as- 
ymptotic jet width h in terms of /3g and X. The constants yg, and y^ can be de- 
termined from equation (28) once the /3's are found. 

Equations of zeroth-order sli p lines in intermediate T -plane. - To obtain the first- 
order solution it is necessary to first find the equations of the zeroth-order slip lines in 
the T -plane, and To this end, note that figures 3 and 5 show that 

jW Wq(T) = 0 whenever T C 


and 


Wq(T) = -1 whenever Tcy^,^^ 

Upon introducing polar coordinates into the T -plane 

T = pe^^ (35) 


equation (19) shows that 


Wq(T) = i 


p cos 2e - yi ro - ya, 

In — 2- (p cos e + To) 

L ^3 73X5 


1 2 ■>'3 ~ ■>'fi 

e - 7T p sin 20 p sin 0 




>^ 3>^6 


0 < e < 7T 


Setting the imaginary part of this expression to zero and to -1 we find that 


- e) + i p^ sin 20 + (yg - yg) p sin 0 = 0 for T C 
2 


(36) 
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and 


-y^rQ0 +-P^ sin 20 + (yg - yg) p sin 0 = 0 for (37) 

2 

Each of these equations has three solutions. The solution e=Tf, 0^p<°° of 
equation (36) and the solution 0 = 0, 0 ^ p ^ of equation (37) are of no interest to us 
since they correspond to the solid boundaries which are mapped into the real axis in the 
T -plane. The remaining solutions to equations (36) and (37) are, respectively. 


70 - 70 ± Vt^6 ■ '^3)^ ^ 0j^ < 0 < 7T for 


2 cos 0 


P = 






4ygy60 cot 0 


2 cos 0 


0 < 0 02 for TC^q 


y<2' 


whore (since p cannot be imaginary and, hence, the quantities under the square root 
signs must be positive) 0^ and 02 are the solutions of the transcendental equations 


tan 0j^^y0 - 73^^ + 47370(61 - 7t) = 0 for 0j > 0" 

>■ 

9 

tan 02^70 - 73^ + 4737062 = 0 for 02 > 0 ^ 


(38) 


and we have used the trigonometric identity sin 20 = 2 sin 0 cos 0. The proper choice of 
the + or - sign in front of the square root will be discussed subsequently. Hence, in 
view of equation (35) we can write 


T = I yg - yg ± V(70 - 73)^ + 4y0yg(0 - 7r)cot 0 (1 + i tan 0) 


with 


0^ < 0 < 7T for Tcy^ 


( 1 ) 


(39) 
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( 40 ) 


T = - 


>^6 - >^3 ± + 473700 cot ej(l + i tan 0) 


with O< 0<02 for 

When the proper signs are chosen, equations (39) and (40) are parametric equations for 
the zeroth-order slip lines with parameter 0. It is convenient to introduce a new param- 
etric variable a;; 


p for 
p - 77 for 


Then in view of equation (38), equations (39) and (40) can be written as 



T = T*(o)) with -co^ < o) < 0 

for 


(41) 


T = T*(o)) with 0 < o) < CO 

m 

for 


(42) 

where we have put 




T=^(co) = i 
2 

re - ^3 ^ " '^ 3 )^ ^ 

for -co^<co<co^ 

(43) 

and is the solution of the transcendental equation 





(76 - rsf tan co^ + 47376^ni = ® 

for 

^m-0 

(44) 


Note that the solutions of this equation will always be such that > 77 / 2 . The 
corresponds to the location along a slip line where d0/dp = 0 (see fig. 8). In order to 
completely determine the parametric equations of and we must determine 

when the + and - signs hold in equations (41) and (42). First note that 

T-"(0) = 7g (45) 


and 


T-(0) = -73 


(46) 
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In order to proceed further we must consider the two cases yg > and yg < yg 
separately. First consider the case where 7q> This corresponds to the situation 
depicted in figure 5. Note that the conditions (45) and (46) show (see fig. 5) that we must 
have 


T = T+(w) for Tcy; 


( 2 ) 

0 


T = T"(co) for T Cy 


( 1 ) 

0 


when CO is near zero. Now for yg > y^ we will have 


T'*'(co)"-«> as — 

2 


Hence, T'*'(co) will transverse the entire contour in a counterclockwise direction as 
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u) varies from zero to n/2. On the other hand, T“(w) will remain finite when w = -tt/ 2, 
and the curve will cross the imaginary axes when this occurs. In fact, it will re- 
main finite even when o) = at which point T'*' and T" become equal. In order to 

obtain the remaining portion of we must take 


T = T‘^(o)) with < o) < - I 

m 2 

for this will continue 9^^^ smoothly from to infinity. 

Hence, we have 


T”(o)) with < w < 0 for [t| < p. 


m 


T = 


T‘''(a)) with -CO < co < - - for |t( > p 


for 


m. 


T = T'^(co) with 0 < CO < I for T C 9^^ 

2 ^ 


> for yg > yg (47) 


where we have put 


Pm" 


ra - ro 




1 + tan w 


m 


Note that the curve is traversed in a counterclockwise direction when co is increas- 
ing and that the ctirve 9^^^ is traversed in a counterclockwise direction with co increas- 
ing for |t| > p^. 

Precisely the same type of reasonii^ now shows that 


T = T‘(co) with - 1 < CO < 0 for 

2 ^ 


T =• 


T'^(co) with 0 < CO < CO for T < p„ 

Xii 111 j 


[t“(co) with — < CO < CO for ItI ^ p 
2 ™ 


for 


y for yg<yg 


mj 


(48) 


Relation between points in physical and intermediate T -planes. - It follows from the 
first equation (10) that the points in the physical plane (fig. 2(b)) are related to the points 
in the T -plane by 
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f 



dWr 


Cq(T) dT 


dT + constant 


(49) 


Substituting equations (18) and (23) in this formula and using the fact, indicated in fig- 
ure 2(b), that the origin of the coordinate system in the physical plane is to be at the 
point 2 result in 


T 

/T + 1\ ^ (T + yg)(T - yg) 

VT + rJ T 


Formulation of First-Order Problem In Physical Plane 

The mapping T~z defined by equation (50) transforms the upper half T -plane ap- 
proximately into the region of flow in the physical plane. The domain is mapped 
into the shaded region of the physical plane shown in figure 9. The curves and 

are mapped into the dashed boundaries and respectively, of this region. 

This region, of course, differs from the true interior of the jet whose boundaries 
and are indicated by the solid (ciirved) lines in figure 9. 
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Relation between quantities on unknown slip lines to quantities on known zeroth-order 
slip lines. - Now the first two groups of boundary conditions in equation (8) are specified 
on the curves and in the physical plane, whose shapes are not known at this 
stage of the solution. As explained in reference 1, however, these boundary conditions, 
correct to terms of order e, can be transferred to and by relating the values 

T TT T TT ^ 


of 


w 


and w^ at an arbitrary point 


of S 

(1) 


(’!') ^ (2) 
^ ' or an arbitrary point z' 


of 


to their values at some neighboring point z^^^ of or z^^^ of respectively, 

by performing a Taylor series expansion of these quantities about z|)^^ or z^q\ Thus, 


(z^j^ - zg^) + . . . for I. = I,H; j = 1,2 

z=zg^ 

- zg^) + . . . for V = I, II; j = 1,2 

where the second and third equations (10) have been used. 

Substituting the asymptotic expansions given by equation (9) into these Taylor series 

3 

and neglecting terms of 0(e ) yields 




>i(4^^) 


\ dz^ 


.(j) 


z=z 


(i) 


+ 0(e^) for = I,II; j = 1,2 


(51) 


w 


(z^J^) = WQ(zg^) + e(w[(zg^) + ?o(4^)^l^ 


+ 0(e^) for 1 ^ = 1, II; j = 1,2 


(52) 


These expressions relate the values of the dependent variables w^, w^\ and at 
the points of the unknown boundaries and to their values on the known bound- 
aries and Sq^\ respectively, with an error of order e^. It can now be shown from 

equation (51) that 

. n \|2 




+ 2 € 


iro 


/?e 


r 

n' 

(4'’) 

5oi 

So 

^ts 3 ^ 




^ol^o 


dz 


,(j) 


z-Zq 


+ 0(e'^) for i^ = I,II; i = 1,2 (53) 
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Boundary conditions alo ng slip lines . - Substituting the expansions (17), (52), 
and (53) into the first two groups of boundary conditions in equation (8) and equating to 
zero the coefficients of e to the first power yields the following first-order boundary 
conditions along the slip lines: 



Thus, equations (54) to (56) are the boundary conditions for the first-order solutions 
on the boundaries S^^^ and ’’transferred" to the zeroth -order boundaries and 

respectively. Hence, the first-order boundary value problem has been transferred 
from one in which the shape of the boundaries is unknown to one in which it is known. 
Note, however, that these boundary conditions involve the variables w ^ , wF ^F, 

(^\ (n\ T TT iili 

z^ , and z'j But and are completely determined in terms of 
w^ and w^^, respectively, by the second two equations (10). In view of this, equa- 
tions (54) to (56) may be thought of as two boundary conditions connecting the variable 

with the variable cF across and (or equivalently the variable w^ with the 
TT u u /T\ T 

variable w ) plus two equations which determine z^^ and z^^ once is known. 

Thus, it follows from equations (55) and (56) 


jW j = 1,2 


(57) 


Hence, in view of the second two equations (10), equations (54) and (57) are boundary 
conditions on and which connect the solution F in with the solution fF 

Ti t; u 1 (2) i- 

in D . And equations (55) and (56) with v = l serve to determine z\ and z\ ' 

I Cl) C2) ^ ^ 

respectively, once is known. (Actually, z^ and z)^ will be determined in a 

slightly different fashion. ) The conditions (57) can be differentiated with respect to Wg 

along the zeroth-order slip lines to obtain, upon using the second and third equations (10), 





for j = 1,2 
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Multiplying this by i and adding it to equation (54) yields the following complex jump 
conditions along and 





0^6 



for j = 1,2 


( 58 ) 


Boundary conditions along solid boundaries. - The bovmdary conditions for the re- 
maining (solid) boimdaries are easily deduced by substituting the first two asymptotic ex- 
pansions in eqviation (9) into the boundary conditions (7) and (8) to obtain 


9*k ?i(z) = 0 for z C 45 and z C 56' 

9*k ?^(z) = 0 for z C 12 and z C 67^ 


arg ^ j(z) = 7 tX for z C 




J1 


arg (z) = ttX for zC 


34 


23. 


(59) 


(z) 0 for z — point 1 and z — point 7 

The second boundary condition (59) can be combined with the zeroth-order boundary 
condition (16) to obtain 


(60) 


9m 


K\{7.) 

?n(z) 


= 0 for z C 34 


9m 


= 0 for z C 23 


Similar remarks apply for the first boundary conditions (59) and (14). Thus, the 
boundary conditions (59) can be rewritten as 


9m 


d(z) 




= 0 for z C 34, z C 45, and z C 56 


(61) 
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= 0 for z C 12, zC23, z C 67 


(62) 


and equations (17) and (60) can be combined to yield 

??(z) 

—0 for z— point 1 and z —point 7 (63) 


Solution of First-Order Boundary Value Problem 


Boundary v alue problem in T -plane . - The boundary conditions (57) and (61) to (63) 
completely determine a boundary value problem (or more precisely, two boundary value 
problems connected along the ciirves and S^^^) for a holomorphic function in the 
zeroth-order region of flow in the physical plane. Under the change of variable z — T 
defined by equation (50), this boundary value problem can be transformed into one in the 
upper half T-plane (fig. 5). The boundary conditions in the T-plane become 


?\(T) 


??(T) 


?0(T) 


^for T 
Jot 


9m. 


?5(? + iO) 

-i = 0 

+ iO) 


for -yg < ^ < yg 




??(? + iO) 

9m =0 for ^ < -y„ and ^ ^ yo 

+ iO) 


??(T) 

— 0 for T ^ 


(64) 


Use of sectio nally analytic fimction 0. - Clearly, the domains of definition of 

and are and respectively. It is convenient to work with the sectionally an- 

T TT 

alytic function © defined in the upper half T-plane in terms of and by 
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0(T) = 


' I 1 I 

©\t) = -= for 

?q(T) 2 


0^^(T) = — 


??(T) U 

-i for TC 9^ 


UT) 


It follows from the boundary conditions given in equation (64) that 0 must satisfy 
the boundary conditions 


0l(T) - 0°(T) = ii(T) 


for TC^X 


for 


0(^ + iO) = 0 for -oo < ^ < 


0(T) - 0 for T 


where we have put 


m(t) = - 


|?0(T)|‘ 


Since ^ ^ and can be no more sir^ular than if the asymptotic expansion is to be 
tiniformly valid, it follows that 0 must be boimded in the entire upper half T -plane in- 
cluding the real axis. 

Now the Schwarz reflection principle shows that the boundary condition (67) will be 
satisfied if, and only if, the function 0 is analytically continued into the lower half 
plane by the formula 

0(T) = 0(T) for ?7 < 0 (70) 


This can be done by first extending the curves and into the lower half plane in 
such a way that they are symmetric with respect to the ^-axis as shown in figure 8. 
This extension can be accomplished algebraically simply by extending the range of the 
parametric variable co so that it is symmetric with respect to the origin in the param- 
etric equations (47) and (48) for the curves and . Thus, the parametric equa- 
tions of the extended curves and shown in figure 8 are 
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T 


and 


rT“(w) with -Wjjj < o) < o)j^ for (t| < 

(t'^M with -Wjjj <w<-I and for It] > 

T = T'''(w) with - 1 < w < - for T 

22 “ 


for Tcy 


( 1 ) 

0 


(71) 


for ye>y3 

i.^1) 


T = ■ 


T = T"(o)) with - 1< O) < - for 
2 2 

|t+(w) with s o> £ for 1 t| < p^ 

1t‘(cd) with -oijjj s o) < - - and 1 < co s for |t| a 

2 2 


for rc^\ 


,(2)f ^6 <^3 


0 


The arrows in figure 8 indicate the direction in which the curves are traversed with co 
increasing. 

Next it is necessary to extend the definition of the jump /u(T) to the negative half 
T -plane. It is clear that this jump must have the same symmetry properties as 0(T) 
with respect to the real axis. Hence, we must have 


p(T) = p(T) for j] :< 0 

Substituting equation (23) into equation (69) shows that for 77 0 

/ 't!^ + 2/?eT _ 1 

But since 


m(T) = i 
2 


(73) 


(74) 


/ !t 1^ + 2/50T + i_Y _ j ^ / !t|^ + 2^e;r^Hj_Y _ J 
\It| + 2>'^^e T + y|/ \ 1 T I ^ + 2y^ /?eT + y^/ 

we see that the symmetry conditions (73) will be satisfied if we take equation (74) as the 
definition of q(T) in the negative half plane also. We now require that the jump condi- 
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tion (66) hold along the extended curves and 9'^^ with fi(T) defined in the entire 
T -plane by (74) and that equation (68) hold in the negative half plane also. Then the 
boundary condition (67) will be automatically satisfied (see ref. 4). 

Determination of 0 by Cauchy integrals alor^ slip lines . - Since /i(T) —0 as 
T the sectionally analytic fvmction 0 which satisfies the jump conditions (66) on 
the extended curves and which has the behavior given by equation (68) at 

T = “, and which is analytic at every other point of the T -plane can be represented by 
two Cauchy integrals as (ref. 4) 



where the direction of integration is such that the region I is to the left when traversing 
the contour. In order to evaluate these integrals we must first transform them to in- 
tegrals with respect to the parametric variable w. 

Since the variable of integration t traverses the contours and this 

change of variable is accomplished simply by substituting in the parametric representa- 
tions (71) and (72) of the contours and When this is done and the directions of 

integration indicated in figure 8 are taken into account, we obtain 


0(T) 


1 

27Ti 



M+(g;) dT+(g;) 

T‘^(o;) - T dc4) 


do) + — 
27ri 



M~(g)) 
T"(w) - T 


dT~(g)) 

dw 


dco 


(76) 


where we have defined M*(w) by 



for 


1 <7 


M"^(w) = < 



2 


[-m(T-"(co)) 

for 

7T 

2 

< 1^1 ^ ^m| 

^ for yg>r 3 

M"(o)) = m(T"(w)) 

for 

1 

S 

3 

VI 



and 
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and 




M^(cu) = fi(T‘^{oo)) for |o;| < co. 


m 




|ju(T"(co)) for |a)| <1 

2 

U|u(T"(a))) for 1< |o)| 


f ^6<^3 


where the minus signs take into account the reverse direction of traversal on certain 
portions of and These formulas can now be combined into the single formula 


= p(T*(co))sgn -± |ct)| sgn(y 3 - yg) 
-2 


for -CJ < CO ^ Ct) 
m m 


(77) 


where 


sgn X = +1 X > 0 

sgn X = -1 X < 0 

In view of the mapping T — z defined by equation (50), this completes the solution of the 
problem since equation (65) determines and in terms of the known fimction 0. 

Expression of first-orde r potential in terms of 0 . - It is also convenient to have an 
expression for w^ in terms of 0. In order to obtain this , note that it follows from the 
first and second equations (10) that 

^I ^ dwj ^ dwj ^ 

Si Sn 

^ dz dT dwg ^ 


Hence, 


dw^ ^ dwg 
dT dT 


and equation (65) shows that 


dwj 

dT 


0^T) —2. 

dT 


2 dT 


( 78 ) 
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Note that z = when T = and therefore equation (6) implies that both Wj(T) 

and Wq(T) vanish when T = -yg. Hence, integrating equation (78) between T = -yg and 
T yields 


Wj(T) = 



©Vt) ^ 


dT + i Wn(T) 
dT 2 ^ 


(79) 


Volume flow in jet . - Let Q denote the dimensionless volume flow per unit width 

through the jet. It follows from the definition of the stream function that, to within an 
2 

error of 0(e ), 




Wo(-yg + iO) + ewj(-yg + iO) 




m 


Wo(ye + 


On using equation (19) for WQ(yg + iO) and the fact that 
T = -yg, this becomes 


iO) + ew^(yg + iO) 


+ O(e^) 


w^(T) and Wq(T) both vanish at 


Q = 1 - e Wj(yg + iO) + O(e^) 


Substituting equation (79) into this expression and using equation (19) again shows that 


Q = 1 + - - e 
2 



©V) — ^dT + O(e^) 
dT 


On substituting in equation (18), this becomes 


Q = 



T (T + yo)(T - yg) 2 

0\t) ^ ^dT + O(e^) 


(80) 


In view of the singularity in the denominator, this integral must first be carried out over 
the path of integration shown in figure 10 and then the limit £q 0 can be taken. Per- 
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V 


3 

-^3 



^6 


Figure 10. - Path of integration for equation (80) in T-plane. 




forming these operations yields 


1 

^ys'>'6 



. (T + r^)(T - rg) 
0^(T) — — — dT 


1 



T (4 + - rc) I, ^ 

0^(^) e L + i0\o) 


Since 0^(4 + iO) is real, equation (80) becomes 


Q = 1 + i + e0^O) + O(e^) (81) 

2 

Since the flow far downstream of the jet becomes uniform and since lim = 1, it 

follows from equation (2) and figure 2 that lim C (z) = wl + Therefore, continuity re- 

X-*oo » 

quirements dictate that the dimensionless net volume flow through the jet Q be equal to 
6yi + €. Upon expandir^ this in powers of e we find with the aid of the expansions in 
equation (12) that 


Q — (l + e6j+. 


. )(1 + 1 e + 
2 


or 


Q = 1 + e 



(82) 
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Hence, equating like powers of e in equations (81) and (82) shows that 


= 0 ^( 0 ) 


( 83 ) 


Since equation (50) sets up a one-to-one correspondence between points in the phys- 
ical plane and points in the T-plane, it is clear that equations (65), (75), and (79) can be 
used to compute the first-order perturbation to the velocity and stream function at each 
point of the physical plane within the jet. 

Because of the fact that once the shape of the jet is known it is quite easy to sketch in 
the streamline patterns, the most important quantities to be obtained from the analysis 
are the shapes of the curves and (see fig. 2) in the physical plane. 

Coordinates of slip lines. - In view of the one-to-one nature of mapping involved, if 
= 0, then z must be a point on the streamline that passes through the point 3 in 
figure 2. In addition since the velocity potential is increasing in the direction 3—1 
along the slip line it is clear that, if /?e w^(z) > l^e w^(Ze"^’’^^) = 0, then z must 

be a point on Now it is easy to see from the zeroth-order solution that 

jWwg = 0 and Wq > 0 for any point These considerations 

imply that the point z^^^ = + ez^^^ will be on the first-order position of if we 

let z^^^ satisfy the equation 


4'"] -o[4'>] ^ I) 


(84) 


This equation shows that when z^^^ = w^[z^^^] = 0 and that Re w^jz^^^j 




Hence, the point z'^' traverses the slip line (to first order in e) as z 


«> as 


traverses the zeroth-order slip line S 


( 1 ) 

0 


,(1) 

0 


Analogous considerations show that the point z^^^ = Zq^^ + ez 
by the equation 


which is determined 




I, 

w (-r 


ih) 


= St 


(2)- 


Wp(-r - ih) 


1 + i 


(85) 


42) 


traverses the slip line to within first-order terms as the point zl^^ traverses the 

( 2 ) ^ 
zeroth-order slip line '. By substituting equations (84) and (85) into the expansion 

given in equation (52), the first-order distances z^^^ and z^^^ from the zeroth-order slip 

lines to the actual slip lines are found to be, respectively. 


41) - 


Cn(T) 


w'(T) -iw(,(T)] 


for TC.^, 


( 1 ) 


(86) 
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z(2) = _ p^(T) - w\(rg + iO) - iwg(T) + i WQ(rg + iO)J for T (87 

where the fact has been used that the curves and in the physical plane are the 

conformal images under the mapping T -*z defined by equation (50) of the curves 
and respectively, in the T -plane. This also shows that 

= z(T) + + O(e^) for (88 

z^ 2 ^ = z(T) + ezf\T) + O(e^) for T CyJ^) (89 

Substituting equations (86) and (87) into equations (88) and (89), respectively, and 
substituting equation (79) into the resulting expressions yield 


z^^^ = z(T) 


L 


Mi 

(T) — y dT for T 

dT ^ 


z^^^ = z(T) 1- 

?n(T) 


0^(T) ^ dT for T ^ 
dT ^ 


where the contour of integration can be taken along and provided that 0^(T) is 
interpreted as the limiting value of 0(T) as T approaches and yj^^ from within D^. 

Static pressvme distributions in jet and outer flow. - It follows from the first two ex- 
pansions in equation (9) that 

C^(z) = ? (z) + eC[(z) +0(e^) = Co(z^| + 2e /?e ^ + 0(e ) for i^ = I,II 

Substituting equation (65) in this equation shows that the magnitude of the velocity is given 
2 

to terms of order e by 

?^(T) ^ t ^ 2e0^T)] + O(e^) for T C (92) 


^^(T) ^ = ^g(T) ^[1 + 2e0^(T)] + O(e^) for T C: 
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It follows from Bernoulli's equation and equation (1) that within the jet 


P - Poo 



= 1 + e - 



2 


and within the main stream 


P - P 


oo 



= 1 - 



2 


Hence, substituting equations (92) and (93) in these equations gives the following expres- 
sions for the static presstire distributions; 


P - Poo 

1 .j2 


= (1 + e) 


1 - 


)(T)|‘ 


- 2e 


Cq(T) ^ 0^(T) for T C 


P - Poo 

1 ,r2 


= 1 - kn(T) 


- 2e 


?n(T) 


^ 0^^(T) for T C 


(94) 


(95) 


All necessary results have now been obtained. However, it is useful to rewrite some 
of these in more explicit form. 


Explicit Formulas for Calculating Quantities of Physical Interest Along Boundaries 

We shall now find the limiting values of 0(T) in equation (76) as T approaches the 
contours and In order to do this we apply the Plemelj formulas (ref. 4), be- 

ing sure to take into account the reversals in the directions of integration indicated in 
figure 8. Once this is done and the parametric equations (47) and (48) of and 

are inserted into the results, we obtain 
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0(T) = 


0"(co) with o) < 0 for |T| ~ 


0’^(o)) with -CO < CO < - — for |t| ^ p 


©(T) = ©■^(w) with 0 < 0 ) < I for T€ 5^1^^ 


for 


for yg>r 3 


©(T) = ©"(o)) with - 1< w < 0 for TC9’iP 

2 0 


©(T) = 


V(w) 

with 

0 < OJ 

for 

|T| 

-Pi 

0'(w) 

with 


for 

|t[ 

-Pi 


for Tcyj, 


for Vc < 




where we have put 


0^{lo) = i /i[T''=(a))] + -L 


T"f(S) - T^(w) dw 


T^(S) - T*(oj) dw 


The difference in sign between the jj.{T^) term and the is due to the reversal in 

the direction in which the contour is traversed with increasing co in certain regions. 

The fact that the proper sign is used in the (l/2)p term can be verified by taking the 
limiting values directly in equation (75). 

Upon differentiating equation (43) we find that 


dT*(u)) i ’"e " " '^ 3 )^ ^ ^ ^'^6^3 ^ 


do) 2 , 

cos (V 


2 cos^ o) sin^ co ^{yq - ^ 


Coo rdinates of slip lines . - Now substitutir^ equations (18), (50), and (96) into equa- 
tions (90) and (91) and using equations (47) and (48) and figure 2(b), we find the following 
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parametric equations for the slip lines 9 "^^^ and correct to first order in e 




+ z"(w,0) + — — z“(o),0) for < ct) < 0 

0 ?-(o>) ” 




>for re>y3 


= -(r + iW + Z^(OJ,0) + — ^ z'^(w,0) for 0 < co < i 


?5(o>) 


lO-; 


z^^^ - le + z;(o), 0) + — - — z (a.', 0) for - - < a' < 0 
° 2 


,( 2 ) 


-(r + ih) + zt(ai, 0) + — - — z'*'(ai, 0) for 0 < o) s cd 

° Co^a-) ” 


[-(r +ih) + z$(a)„,0) + z-(o;,a'„^) + ^ V] 

where we have put 


for - < 0 ) ^ o)„ 
2 


^ for yg < yg 


z±(c^l,Wo) = - 





L- " >-3irT*(“) - r el mfM 

T*(w) 




do) 


dci) 


0=^(0)) 


[t^M . r3][T"(c.) - ^ ^ 


T*(w) 


do) 


and 


?S(co) = 


T="(cc) + 


nx 


T=^(w) + 1 


Summary of Equations 


For convenience, the most important equations are now summarized. 


( 100 ) 


( 101 ) 


(102a) 


(102b) 


( 103 ) 
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Relation between )3-values and dimensionless length of deflecting plate; 


I = 




Jr(^3 + %)(^e " 



j (/3g - (r)(/3g + (j) 
(7^ (l + a)^(/35 + a) 


do 


( 31 ) 


Relation between /3-values and dimensionless orifice offset distance: 


r = 


(1 - hY 


77(^3 + %)(^e - ^5) 





I + o){a - /3g) 
0^(0 -/3g)(l -cr)^ 


da 


(^5 " 


1 - — - - 1) 


, ^ . J) . _ ,y 


(33) 


(32) 


Dimensionless orifice width (reciprocal of dimensionless asymptotic zeroth-order jet 
width): 

h = J- 




(34) 




^^3 


^6 


1+^3 
^4 = ^5 
^6 ■ ^5 


1-^c 


J 


(yg - ^3)^ tan + 4y3rgO)^ = 0 for o;^ - 0 


(28) 


(44) 


T"^(a)) = -jyg - V3 ± V(^6 ~ ^3)^ w (1 + i tan w) for -co^ < co < o)^ (43) 
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(99) 


dT"‘'(ct)) i ^6 " ’^3 ^ ''■ ^(sin co cos coe^'" - coe"^'") 


,11^ ,.„-lCt^\ 


dcc) 2 2 

2 cos" w sin" w 


,2 .. _,_2 


Vt^6 - ^ 3 )^ + ^ygrgw cot 


w 


m(T=") - i 
2 


iT^I^ + 2 ygeT^ +. 1 

[\Jt^|2 + 274 ke T* + y\^ 


- 1 


M'*'(a)) 


ju[T''=(c<;)] sgnj^^i jw| sgn(yg - yg) 


Q^{u>) = i /i[T'*'(a))] + 


Cl) 

/» m 

^ / 

jTTi / 

«/_r.) 


M^(w) dT^(S) 


27Ti / T’^(o)) - T=^(w) dco 

-Cl) 

m 


27Ti 


I 


in 


M^(w) dT*(w) 


T"^(a)) - T^(a)) dw 


dco 


0(^ + iO) = Re — 


771 


/ 


Cl) i^CL) 

m . / m 

+/ 


M~^(g)) dT'^(g^) ^ j M~{w) dT~(gj) 

T"^(c 4)) - ^ dco I T“(w) - ^ do) 

L*'0 *'0 


z±(a)i,o)Q) = 





T^(w) + yg^ T=^(c4j) - yg] 


(co) 


T*(w) 


do) 


dcu 


z*(WpO)o) 





Q^iu>) 


[t*(") * - 5-e] dT= 


'{oj) 


T*(w) 


dco 


dco 


r“rpi 


?q(") = 


T^(co) + y4 
.T*(w) + 1_ 


(74) 


(77) 


(98) 


(108) 


(102a) 


(102b) 


(103) 
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Coordinates of the slip lines: 


,(1) = 


.(2) J 


= + Zn(o),0) + — 5_z"(c(j, 0) for - _ < u> < 0 

° ^ 

l(r + ih) +zi(a),0) + z^{u,0) for 0 < to < £*7^^^ 

CqM 

-(r + ih) + Zq{‘^^, 0) + Zq^co, coJ + ^ ^ for 2 < to > 

^0 ^ ^ 


( 100 ) 


>• for Yd < 


’>'6 ^^3 


7.(1) = 


+ zr(to,0) + — i — z'(to,0) for -to s to < 0 
" ?-(to) “ 

f 4- ZQ(-to^, 0) + zJ(to, -to^) + ^ 0) + z+(to, -to^)] 

Ca(w) 


( 101 ) 


for -to„-'toV.^ 
■" 2 


Wor 


7(2) 


-(r + ih) + zi(to,0) + -— -- z*(u.’,0) for 0 • w •' - 
° Co(to) 2 


OUTLINE OF CALCULATIONAL PROCEDURE 

The numerical evaluation of the analytical results was carried out by using complex 
arithmetic. Hence, there was no need to analytically separate the real and imaginary 
parts of the various equations which were used. In order to evaluate the final expres- 
sions, certain parameters and results which are completely determined by the zeroth- 
order solution must be found. The procedures involved in accomplishing this will there- 
fore be discussed first. 

As shown by equations (28) and (31) to (34) (refer to the section Summary of Equa- 
tions), the physical quantities I, r, and h are related implicitly to the parameters (3^, 
/3g, and /3g in a complicated manner; and the parameters y^, y^, and yg are in turn 
related to the /3's. For a given deflection angle A, = a/v, it is convenient to choose 
values for the jet contraction ratio Aq/H = l/h and for /3g. Then from equation (34) the 
parameter /3g is fomid, and the parameters y^ and y^ are obtained from the first two 
equations (28). Equation (32) can then be solved for /3g; and with all the j3's known, Z, 
r, and yg are obtained from equations (31), (33), and the last of equations (28). Since 
l/h and convenient computer input parameters, grids can be prepared as shown 

in figure 11 to aid in choosing input parameters from the physical data. Note that the 
magnitude of the orifice width H is arbitrary. The values in figure 11 are for positive 
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(a) Plate angle o, -3CP. (b) plate angle a, -60°. 

Figure 12. - Relation of jet contraction ratio and mapping parameter ^3 to geometric quantities for negative values of angle a. 


values of the angle of. Some results will also be given for negative a and zero orifice 
offset distance (R = 0), with the relations between Aq/H, and L/H for these cases 
given in figure 12. 

Now that the physical parameters have been expressed in terms of the y's, the co- 
ordinates of the jet bounding slip lines can be obtained. The transcendental equation (44) 
is first solved for co^. With the y's known, the functions T*(oi) and dT^/du) can be 
evaluated from equations (43) and (99) and then iU(T^) and M*(co) can be obtained from 
equations (74) and (77). These quantities are used in equations (98) and (103) to obtain 
Q^{w) and ?q(w). Then the integrations in equations (102a) and (102b) are carried 
out to yield the functions Zq(o 1 j^, oiq) and z*(co^, Wq). For various values of the total 
pressure difference parameter e, the coordinates along the slip lines can be found from 
equations (100) and (101) depending on whether y^ is larger or smaller than yg. 

DISCUSSION 

The objective of this report was to develop a technique to analyze the flow interaction 
when a two-dimensional jet is injected between two moving streams. Such an analysis 
must account for conditions where the total pressure in the jet differs from that in the 
outer stream. The analysis must also account for the variation of static pressure along 
the jet bounding slip lines, a feature not taken into accoimt in classical free-jet analyses. 
The analysis was applied to the configimations shown in figure 1, and the figures that are 
presented in the following discussion will show typical slip line patterns. In all the fig- 
ures the geometry has been nondimensionahzed by the orifice width H shown in figure 2(a) 
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Positive Deflection Angles 


First, consider the results for positive values of the angle a. Two values, a = 30° 
and 60°, were considered; and for each of these, results are presented for three dimen- 
sionless plate lengths L/H = 1, 2, <and 3 and two orifice offset distances R/H = 1 and 2. 
The results are shown in figures 13 and 14, which correspond to a = 30° and 60°, re- 
spectively. Table I gives the specific values of the parameters Aq/H and /Sg for these 
cases. 
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X 
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(c) L/H = 3; R/H = 1. 

Figure 13. - Jet slip line boundaries for outlet with 30P plate deflection angle. 
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TABLE I. - ZEROTH -ORDER CONTRACTION RATIO AND /3g AS 


FUNCTION OF GEOMETRIC PARAMETERS 


Plate 

Dimensionless 

Dimensionless 

Zeroth-order 

Mapping 

Figure 

angle , 

plate length, 

offset 

contraction 

parameter, 


a 

L/H 

distance, 

ratio. 

h 





R/H 

Aq/H 



30 

1 

1 

0. 67673 

0. 63568 

13(a) 



2 

1 

. 59198 

. 64244 

13(b) 



3 

1 

. 54410 

. 65225 

13(c) 



1 

2 

. 75562 

. 67882 

13(d) 



2 

2 

. 68060 

. 67645 

13(e) 



3 

2 

. 63411 

. 67821 

13(f) 

60 

1 

1 

.40540 

.42008 

14(a) 



2 

1 

. 27970 

.42814 

14(b) 



3 

1 

.22175 

.44116 

14(c) 



1 

2 

. 53905 

.46570 

14(d) 



2 

2 

.41110 

.46160 

14(e) 



3 

2 

. 34041 

.46369 

14(f) 

-30 

1 

0 

1. 57478 

1. 14683 

15(a) 

-30 

2 

0 

1. 76227 

1. 18917 

15(b) 

-30 

3 

0 

1. 88322 

1. 22134 

15(c) 

-60 

1 

0 

2. 17472 

1. 70739 

16(a) 

-60 

2 

0 

2. 60760 

1. 74679 

16(b) 

-60 

3 

0 

2.90804 

1. 77788 

16(c) 


The ratio Aq/H is the zeroth-order contraction ratio of the jet. The results in 
either table 1 or figure 11 show how this contraction ratio decreases as either L/H is 
increased (with R/H fixed) or R/H is decreased (with L/H fixed). Increasir^ the 
angle a also causes Aq/H to decrease, as would be expected from the geometrical 
configuration. 

Figures 13 and 14 show the slip lines for various values of the total pressure differ- 
ence parameter e. It follows from the nature of the asymptotic expansions in equa- 
tion (9) that the parameter e cannot be so large as to cause a large deviation in the slip 
lines from their zeroth-order positions. As e is increased, the total pressirre in the 
jet is increased above that in the stream, as shown by equation (1). This produces two 
effects on the slip line configirration: the jet width expands as e increases; and the jet 
deflection tends to persist over a lor^er path, that is, it is more difficult for the outer 
stream to turn the jet into the horizontal direction, resulting in an increase in jet 
penetration. 
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results depend on the flow remaining attached to the upper boimdary by means of the 
Coanda effect, and this attachment can be promoted by boundary layer control. 

All the results in figures 15 and 16 are for zero offset of the orifice lip (R/H = 0). 
The contraction ratio Aq/H (see table I) is larger than unity in this instance as the flap 
is expanding the jet by virtue of the attachment of the upper jet streamline. When the 
value of e is increased, the width of the jet tends to expand and its path is more 
horizontal. 


CONCLUDING REMARKS 

An analytical technique was devised for determining the flow characteristics of an 
inviscid two-dimensional jet interacting along both of its boundaries with an outer main- 
stream flow. The velocity variations in the outer flow produce a priori unknown static 
pressure variations along the jet boundaries, a condition not accounted for in classical 
jet analyses, which ordinarily have a uniform static pressure (stationary outer flow) 
along the jet free streamlines. The analysis also accounted for the total pressure in the 
jet differing from the free -stream total pressure, with the restriction that this difference 
does not become too large in comparison with the free-stream kinetic energy head. 

Specific results were carried out for two config\arations of the channel used to inject 
the jet into the stream, and the shapes of the bounding jet slip lines were calculated. 

The results show how the jet penetration and contraction or expansion depend on the total 
pressure in the jet relative to that in the free stream. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 23, 1972, 

136-13. 
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APPENDIX - SYMBOLS 


b characteristic length, put equal to Aq 

Cq integration contour, fig. 7 

D^,D^ flow regions in physical plane 
Dq, Dq zeroth -order flow regions in physical plane 

flow regions in T-plane 
H orifice width 

h dimensionless orifice width, H/Aq 

L length of plate at angle to flow 

I L/Aq 

defined in eq. (77) 

O order symbol 

total pressure in main stream 
Pj total pressure in jet 

Y. Cauchy principal value 

p static pressure 

p^ static pressure far upstream from orifice 

Q volume flow through jet 

R orifice offset distance 


r 

g(l) g(2) upper and lower slip lines (in physical plane) 

zeroth -order position of upper and lower slip lines 
y(l) y(2) gjjp Hues in T-plane 
T intermediate variable, T = ^ + irj 

T* defined in eq. (43) 

U X-component of velocity 

^ U/V^ 

V Y -component of velocity 

free -stream velocity 


Ill 


V 

w 

X 


X 

Y 


y 


,(i) ,(2) 


a. 

B(m,n) 

r 


y3,y4,re 

A 

6 


€ 


v/v„ 

dimensionless complex potential, W = (p + ixl/ 

coordinate in physical plane 

X/A(, 

coordinate in physical plane 

Y/Ao 

dimensionless complex physical coordinate, z = x + iy 
defined by eq. (102) 

dimensionless coordinates of points on slip lines 

angle of plate relative to undisturbed main stream 

beta -function 

parameters in -plane 

gamma-function 

parameters in T -plane 

asymptotic jet width 

A/Aq 

(Pj - P J/(l/2)pv2 


"O 

? 

V 
0 
0* 

e 

01 , 02 
X 

M 

V 

P 


radius for contour integral in figs. 7 and 10 

dimensionless complex conjugate velocity, ^ = u - iv 

coordinate in T -plane 

defined in eq. (65) 

defined by eq. (98) 

polar coordinate in T -plane 

solutions of trancendental eqs. (38) 

a/n 

defined by eq. (74) 
dummy index 
coordinate in T -plane 
density of fluid 
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a coordinate in S2 -plane 

T coordinate in -plane 

# velocity potential 

^ ^/^O^cc 

'Jr stream function 

Q intermediate mapping variable, = a + ir 

ct) parametric variable defined after eq. (40) 

c^m solution of trancendental eq. (44) 

Subscripts: 

0 zeroth-order quantity 

1 first-order quantity 
Superscripts : 

I inside jet 

n in main stream 

1 on S^^^ 

2 on S^^^ 


complex conjugate 
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